ON SELF-SIMILAR COLLAPSE OF DISCONTINUES DATA FOR THIN 
FILM EQUATIONS WITH DOUBLY DEGENERATE MOBILITY 



u(x, 0) = sign a; 



V.A. GALAKTIONOV 

Abstract. As a basic model, the fourth-order quasilinear thin film equation (the TFE- 
4) with a concentration-dependent doubly degenerate mobility coefficient vanishing at 
two equilibrium levels u = ±1, 

u t = — (|1 — u 2 \u xxx ) x , 

is studied. The basic Riemann problem for the TFE-4 with discontinuous initial data 

1, x > 0, 
-1, x < 0, 

is considered. This problem is shown to admit self-similar solutions of the form 

u + {x,t) = f(y), where y = x/ti and / solves the ODE 

-(|l-/ 2 |/'")' + i/'2/-0, /(±oc)=±l. 

The similarity profiles are different for the zero contact angle and zero-flux FBP in a 
bounded domain and for the Cauchy problem in K x M. + . The similarity profiles can be 
also obtained by branching at n = for the mobility coefficient |1 — f 2 \ n that involves 
the linear bi-harmonic equation u t = —u xxxx . 

Finite propagation in the same Riemann problem is studied for the TFE-4 with the 
unstable diffusion term such as 

u t = -(|1 - u 2 \u xxx ) x - (|1 - u 2 \ m u x ) x , 

where m > is a parameter, and dynamic interface equations are derived. Similar 
conclusions on similarity solutions of Riemann's problem apply to the 2mth-order TFEs, 

u t = (-l) m+1 (\l-u 2 \ n Dl m - 1 u) x (m> 2 integer, n > 0). 



1. INTRODUCTION: RlEMANN'S PROBLEMS AND SIMILARITY SOLUTIONS 

We study the collapse of initial discontinuities (singularities) in higher-order quasilinear 
degenerate parabolic equations of thin film type. Such singularity problems are classical 
in PDE theory and systematically occurred in entropy mathematics of conservation laws 
and hyperbolic systems in the 1950s. 

1.1. On Riemann- type problems in nonlinear evolution PDEs. Riemann's prob- 
lems play the truly special role and are key for many nonlinear evolution equations with 
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partial derivatives, to be discussed in a short survey below. Typical Riemann's problems 
have been studied for nonlinear PDEs with quasilinear spatial differential operators of the 
orders up to three that are denoted below as equation classes (1), (2), and (3). 

(1) Conservation laws. This is the most classic matter associated with gas dynamics 
problems, where the shock and blow-up Guderley-type regimes were actively studied since 
the 1930gj. Namely, Riemann's problems are well-known for scalar conservation laws such 
as the first-order ID Euler equation as the key representative, 

(1.1) u t + uu x = in lxl + . 

According to classic theory (see details in Smoller [33], Part III]), for such first (and, odd-) 
-order PDEs, two Riemann's problems are key with discontinuous initial data 

(1.2) S-(x) — ~ sign x and S+(x) = signx. 

The first initial function leads to the shock wave, which is the entropy stationary solution 
u~{x, t) = S-(x) for all t > 0, while the second one gives rise to the rarefaction wave 

— 1 for x < —t, 

1.3) u + (x,t) = f(^) = { f for \x\<t, 

1 for x > t. 



The above formula ( 11. 3D shows the generic similarity mechanism, with the scaling invariant 
j, of the collapse of this non-entropy initial singularity S + (x). For the first data in ( 11. 2D . 
replacing t t— > — t, f i— ► — / in (11.30 . gives the formation of the shock S-(x) in a blow-up 
singularity as t —>■ 0~ via the self-similar solution, with the same profile /, 

(1.4) u-{x,t) = -f{jpz) ->S-(x) as t^O-. 

Riemann's problems describing shock wave formation and collapse of other discontinu- 
ities are also classic for strictly hyperbolic systems since the 1950s; see Bressan [13] and 
Dafermos [17] for history and main results. 

(2) Quasilinear heat equations. The collapse of initial discontinuity S+(x) for the 
second-order parabolic equations such as 

(1.5) u t = (k(u)u x ) x , 

where k(u) > is a given coefficient (k(u) = \u\ corresponds to the signed porous medium 
equation), is also described by similarity patterns with the scaling invariant -t= (sometimes 

called the Bolzman substitutior 



u+(x,t) — f(y), y = where / solves the ODE 
(k(f)fj + | f'y = in R (/(±oo) = ±l). 



1 First results on the formation of shocks were due to Riemann himself in 1858 [40]; see |15) . 
2 Similarity solutions were used by Weierstrass around 1870, and by Bolzman around 1890. In parabolic 
PDEs, this rescaled variable y, together with the backward blow-up one y = — p= for t < 0, was freely 

and effectively used by Sturm in 1836 [IS] in his classification of all possible multiple zeros of u(x,t). 
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To our knowledge, first results on existence of such similarity solutions were due to 
Polubarinova-Kochina (1948) [3H]- Later, such ODEs have been studied in detail for 
arbitrary nonlinearities k(u); see Atkinson-Peletier [HE]- By parabolic interior regular- 
ity, the opposite pattern u-(x, t), describing formation of a shock, is nonexistent, i.e., 
discontinuous solutions cannot appear evolutionary in parabolic PDEs (II. 5p . 

(3) Rarefaction and shock waves for NDEs. Riemann-type problems on discon- 
tinuous solutions S_(x) and collapse of S + (x) are actual for third-order nonlinear disper- 
sion equations (NDEs-3), e.g., such as 

(1-7) u t = (uu x ) xx , 

which, besides its applications in various problems of shallow water theory and compact 
pattern formation (see references in [321 Ch. 4]), are a natural extension of the first-order 
conservation law (11.11) to the third- and other odd-order cases. Actually, writing (11.71) in 
a "pseudo-differential" form 

(1.8) Pu t + uu x = 0, 

with a standard definition of the inverse Laplacian P = (— D^)' 1 > (in a bounded or 
unbounded interval), (11.71) can be viewed as a non-local conservation law. Note that, by 
no means, (11.71) can be treated as a hyperbolic system. 

Concerning the rarefaction wave for (11.71) with data S + (x), this is again described by 
the similarity solution 

(1.9) u+(x,t) = f(y), where y = ^ 

and / solves the following ODE obtained on substitution into (11. 7ft : 

(1.10) iff)" + ify = Q (/(±oo) = ±i). 

This is a more difficult ODE, though general understanding and some results on existence 
of rarefaction (and also shock) profiles are available; see [281 § 7], [321 p. 173], and [29l [30] 
for more advanced <5-entropy theory for the NDEs such as (11.71) and higher-order ones. 

Similarly, finite-time blow-up formation as t — > CP of the shock S-(x) is given by the 
similarity solution 

(1.11) u-(x,t) — f(y), where y = x/(— t) 1 / 3 , 
with the reflected profile —f(y) satisfying the ODE (11.101) . 

1.2. (4) Riemann's problems for thin film models with degenerate mobility. 

Thus, as the next extension of the three classes of PDEs characterized above, we naturally 
arrive at the fourth-order quasilinear parabolic PDEs, for which Riemann's problems, 
though not leading to shock waves, become more difficult. 

In this paper, our main model is the fourth-order thin film equation (the TFE—4) 
with concentration-dependent doubly degenerate mobility coefficient b(u) = 1 — u 2 ; see 
[HI [22l E] for further applied motivations and various results. For convenience, in the 
dimensionless form, we write it as follows: 

(1-12) ut = -(|1 - u 2 \u xxx ) x , 
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so it is degenerate at two levels u = ±1, along which one can expect finite propagation of 
interfaces. We pose ( jl,12p either in a bounded domain (the case of the zero contact angle 
and zero- flux FBP), or in the whole space R x R + (this we call the Cauchy problem, CP). 
Concerning further physical derivations, applied motivations, history, and discussion of 
related models, we refer to Novick-Cohen's most recent monograph [37J. 

The absolute value in the nonlinear coefficient 1 1 — u 2 | in (11.121) is introduced to cover 
the case of the Cauchy problem or other FBPs, where, as we show, for initial data within 
the range, uq(x) G [—1, 1], the solution becomes oscillatory about u = ±1 and gets out 
of the physical range [—1,1]. As usual in TFE problems, this may cause difficulties in 
the physical motivation of such a situation. Nevertheless, such solutions are necessary for 
general TFE theory and for understanding the difference between the Cauchy and free- 
boundary problems for degenerate TFEs. In this case, we are obliged to take the absolute 
value |1— u 2 \ of the mobility coefficients in (I1.12p to keep the equation degenerate parabolic 
for all MGR. Otherwise, the parabolic problem becomes partially backward in time and 
loses its classic local well-posedness. Such instabilities of the mathematical nature can 
be suspicious from the physical point of view, though sometimes should be taken into 
account and lead to Young's measure solutions (not properly developed for such higher- 
order TFEs). 

For the zero-flux and zero contact angle FBP and solutions uG [— 1, 1] from the range, 
one can use the standard TFE-4 



Existence and some regularity results for weak solutions of (11.121) and similar PDEs with 
extra low-order terms are well-known since the pioneering paper by Bernis and Friedman 
[7] ; see also [22] and results and references in more recent papers [TH EE2l [231 H21 HZ] • 

1.3. Two basic Riemann's problems for the TFE— 4 and similarity reduction. 

In view of reflection symmetry x t— > — x of (I1.12p . there exists the single Riemann problem 
for the TFE-4 (fLT2l with initial data 



We call this the Riemann-1 problem or simply the RP-1. 

By parabolic interior regularity, one can expect to have a sufficiently smooth solution 
of the FBP (11.131) . (11.141) for t > as a kind of rarefaction wave in parabolic setting. 
The collapse of such initial singularity and propagation of interfaces for small t > are 
key for understanding multi-phase phenomena within of framework of various nonlinear 
problems of TFE, Allen-Cahn, and Cahn-Hilliard type. 

For the sake of mathematical consistency (and not only that), we also introduce another 
Riemann-type problem with a half of the step in (I1.14p setting, 



(1.13) 



u t = -((1 - u 2 )u xxx ) x . 



(1.14) 




(1.15) 
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where 6(x) is the Heaviside function. We refer to this as the RP-|. Notice the principle 
difference in comparison with the RP-1 in (II . 14p : for fll. 15j) we expect finite propagation 
along the equilibrium level u = 1 only, while on {u = 0} the propagation is infinite. 

Our main goal is to demonstrate that, in general, both Riemann's problems fll . 121) . 
(11.141) and (11.151) admit solutions with the self-similar structure 



1.16) u + (x,t) = f(y), where y 



174 • 



Then for the RP-1, / solves the boundary- value problem for the ODE which, by anti- 
symmetry, we pose in R + only, 

(117) B (/) = -(|l-/ 2 m , + s/ / y = for y>0, 

with two anti-symmetry conditions at the origin, 

(1.18) /(()) = /" (0) = 0. 

More precisely, for the FBP setting, we need to find a finite yo > such that the zero 
contact angle and zero-flux conditions hold, 

(1.19) f(y ) = i, f(yo) = (\i-f 2 \f")(yo) = o (fbp). 

Similarly, for the RP— §, with data (11.151) we obtain the ODE (11.161) in M with the regular 
condition at y = — oo 

(1.20) /(-oo) = 0, 

and the free-boundary conditions (11.191) . 

For the CP, we consider equation (I1.17P with the single condition 

(1.21) /(+oo) = l (CP). 

Actually, we show that the solution of the CP has also finite interface position yo.CP > 
that is larger than that for the above FBP. We then set f(y) = l for y > yo,CP- With this 
continuation, the actual regularity of the CP solutions at the interface point y = yo,CP is 
not known a priori, but the solutions are expected to be smoother than for the FBP, so 
(11.191) are always valid. For instance, we show that, for the TFE-4 with parameter n > 0, 

(1.22) u t = -{\l-u 2 \ n u xxx ) x , 

the actual regularity of CP profiles at the interface is more than C'™' -1 ([•] stands for the 
integer part), so solutions can be arbitrarily smooth for small n > 0, [24]. 

Thus, as a key feature, in the CP, the profile f(y) admits the trivial extension beyond 
the interface: 

f(y) = l for y > y 0iCP . 

For the FBP, such an extension is impossible, and actually, makes no sense, since, by 
definition, the FBP is posed in a bounded moving domain. Therefore, the problems with 
conditions fll . 19[) and (11.211) lead to completely different similarity profiles f(y). 



In general, we expect that, in the FBP setting, the RP-1 fll.17p - fll.19l) has a countable 
set of solutions {/&(?/)}, where fk(y) approaches the unique profile of the CP satisfying 
(11.211) . as k — > oo. The first FBP profile is such that 

(1.23) /i(y)e(0,l) for all y e (0,y ), 

and thus gives the stable (generic) solution of the Riemann problem for the TFE-4. Other 
FBP profiles fk(y) for all k > 2 have larger range, <f_ [0, 1]. Though only the first fi(y) 
profile seems to make much physical sense, we will study others to demonstrate mostly 
mathematical properties of the ODE ( 11. 17ft . 

A similar "FBP-CP interaction" is shown to occur for the RP-|. 

1.4. Connection with the linear bi-harmonic equation for n = 0: towards oscil- 
latory behaviour and branching in the CP. It is curious that structurally, the sim- 
ilarity solution of the CP is close to the corresponding solution for the linear bi-harmonic 
equation 

(1.24) u t = —u xxxx in M x M + . 

In particular, the fundamental solution of (11.241) has the similarity form 

(1.25) b(x,t)=t-$F(y), y = &, 
where F is the unique symmetric solution of the problem 

(1.26) BF = -F^ + \F'y + \F = in R, / F = 1. 

The kernel F = F(\y\) is radial, has exponential decay, oscillates as \y\ — > oo, and, for 
some positive constant D and d — 3- 2~ n / 3 (see [2U P- 46] and more details in Section Hj), 

(1.27) \F(y)\ < De"* |4/3 inf. 

Therefore, by convoluting fll.25p with data (j!.14p . the following solution to the RP-1 is 
obtained: 

, y +oos 

(1.28) u(x,t) = f (y) = b(x,t)*signx=( J - J )F(z)dz, y = ^. 

— oo y 

Hence, as earlier, /(0) = f"(0) = 0, and f(y) is oscillatory about ±1 as y — > ±oo. 

The linear PDE (11.241) is obtained from the more general TFE-4 (11.221) by passing to 
the limit n — > + . Such a continuous (homotopy) connection for some classes of solutions 
yields that the behaviour of solutions near u = ±1, which for n = is oscillatory (there 
exist infinitely many intersections with ±1), is inherited by solutions of the TFE-4 (11.221) 
at least for sufficiently small n > 0. Actually, the oscillatory changing sign behaviour for 
TFEs takes place in the parameter interval [24"l § 7.2] 

(1.29) ne[0,n h ), with n h = 1.758665... , 

where at n = nh there occurs a heteroclinic bifurcation in some related nonlinear ODE. 

The existence and uniqueness of the solution fo(y) of Riemann's problem for n = 
given in (11.281) suggests another approach to the nonlinear one (11.171) . Namely, in Section 



HI we intend to show that there exists a continuous n-branch of such similarity profiles 
for the TFE-4 (11.221) . which is originated at n = from f . 



1.5. On extensions. We must admit that Riemann's FBP and CP for higher-order ODEs 
( 11.171) are quite difficult, and we can give a rigorous justification of existence for the FBP 
only. By a combination of analytic techniques and numerical evidence, we also treated 
the CP and establish that the similarity approach correctly describes the collapse of such 
singularities for more general models. Furthermore, we can extend our similarity analysis 
to 2mth-order TFEs such as (the TFE-2m) 

(1.30) u t = (-l) m+l (\l-u 2 \ n Dl m ~ l u) x (n>0), 

with the discontinuous shock (I1.14p . for which we again pose the corresponding zero 
contact angle, zero-curvature, and zero- flux FBP in a bounded domain and for the 
Cauchy problem in R x R + . In both cases, the similarity solution for the RP-1 is 

x 

(1-31) u+(x,t) = f(y), where y = 

and / solves a boundary-value problem for the ODE which, by anti-symmetry, we pose 
for y > only, 

(1.32) (-l^ll - f\nf(2m-l)y + J-f>y = , /(0) = f"(0) = ... = f 2m -\0) = 0. 

The zero contact angle FBP now involves m + l "zero contact angle" and flux conditions 

(1.33) f(y ) = l, f\y ) = f"(y ) = ... = f^- 1 \y ) = 0, (|1 - f\ n f^)(y ) = 0. 

For the CP, we need just a single condition (11.211) . The case of the FBP is easier 
but still is a hard problem. In the case of the CP for large m's, any geometric-like 
shooting techniques for the ODE (ll.32p are illusive, since we have to deal with truly 
multi-parametric spaces. On the other hand, the n-branching approach at the branching 
point n = (where there exists the unique Riemann's profile fo(y) as in (11.281) ) fits for any 
m and implies existence of a suitable similarity profile, at least, for n > small enough. 

Concerning extensions of the fourth-order TFE (11.121) . we consider the Cahn-Hilliard 
PDE with the extra unstable linear diffusion term 

(1-34) u t = - u 2 \u xxx ) x - u xx . 

However, this PDE is less physically motivated, though presents interesting interface 
properties. Though the collapse for small t > of the initial singularity is described 
by an analogous similarity solution (11.161) . we show that the dynamical laws of interface 
propagation changes essentially in both the FBP and the CP. We derive the dynamic 
interface equation in both cases, which establishes the connection between the interface 
speed s(t) and some interface differential operators. 

Finally, we study the case of nonlinear unstable diffusion 

(1.35) u t = -(|1 -u 2 \u xxx ) x - (|1 -u 2 \ m u x ) x (m>0), 

which emphasizes some interesting mathematical phenomena. 
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2. The FBP: on existence, uniqueness, analyticity, and stability 



Here we consider the RP-1 (07j> - (Ogjl . First, we notice that the FBP and the CP 
have different profiles, and the latter one has faster propagation with larger rescaled 
interface point, 

?/o,fbp < Vo,cp- 

2.1. Local behaviour for the FBP: existence of a unique local smooth solution. 

We need to establish some basic local properties of the solutions near the free-boundary 
points. As we have mentioned, for the FBP we consider solutions in the proper range 
/ G (—1, 1), so that (11.171) takes the form 

(2.1) -((i- AT)' + = o. 

First of all, by classic ODE theory [T6J Ch. I], f(y) is smooth and analytic at any non- 
degeneracy point, at which / ^ ±1. Close to the interface at y = y , the local behaviour 
for solutions of (12.11) is well-known [TOj [26] and is given by the following expansion (still 
formal and hence to be justified): 

(2.2) f(y) = 1 - C(y - yf + § (y - yf + |g (y - yf + O((y - yf), 

where C > is an arbitrary fixed constant. We will show that, actually, (12.21) is the 
expansion of an analytic function. Functions (12.21) compose a 2D bundle with parameters 
{yo, C} to be matched with another 2D bundle at the origin with conditions (jl.l8p . 

(2-3) f(y) = A iy + A 3 y 3 + 0(y 5 ), 

with constants A\ > and A% G R. Actually, the asymptotics (12. 2p give a smooth 
continuous (relative to the presented parameters) penetration of orbits into the space of 
analytic functions f(y) G (—1,1). Moreover, we will show that even (12.21) is the expansion 
of an analytic function at y — yo, so we actually deal with globally analytic orbits. 

Since the papers [TP] |2"6] dealt with third-order ODEs only obtained on integration 
once the divergent fourth-order counterpart (as a hint, see (I2.25P below), we need to 
re-derive the expansion (12.21) for our fourth-order equation (12. ip . to say nothing about 
analyticity, which, to our knowledge, was not addressed in the literature. It is curious 
that a standard application of contraction principle for the FBP then demand analysis 
in functional spaces with singular (unbounded) weights. This approach is rather general 
and applies to 2mth-order ODEs for TFEs as in (11.321) . 

Proposition 2.1. For any C > 0, equation (\2.1\) has a unique local solution for y < yo 
satisfying (\2.2l) . which is strictly monotone increasing on some interval y G (yo — S,yo). 

Proof. We begin by integrating (12.11) over (y,yo), and using the free-boundary flux con- 
dition in (11.191) obtain 

yo 

(2.4) -((i-p)n(y) = \Jfydy. 

y 
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In order to simplify further calculus, let us explain the origin of the unique monotone 
solution of (12 .4p . We represent the right-hand side in the equivalent form 

-((i - p)f "M = \ r; fydy = -\ - fy y \ d y 
(2 ' 5) = -i { - m'dy - -f)M = w-f)y+\ ca - /) **• 

Therefore, dividing by 1 — f 2 = (1 — /)(1 + /) yields the following perturbed equation: 

(2-6) -f" = Kf, y) = h 1 + h 2 = ^ + 4{1 _ f ) (1+f) /;°(1 - /) dy. 

Obviously, the first term is non-singular at the point (y , 1) for / w 1, y « y , so, by 
classic theory, the unperturbed ODE 



v 



( 2 - 7 ) - 4(1+/) 

possesses a unique local analytic solution satisfying, according to the asymptotics (12. 2ft : 

(2.8) /M = l, fM = 0, /" ( W ) = -2C < 0, f"(y ) = < 0. 

Note that any C 3 -solution satisfying (I2.8P is strictly monotone increasing for y m y$. 

The second term in (I2.6P contains certain weak singularity in the same class of smooth 
solutions, so we should treat it in a weighted space. As a standard procedure of Perron- 
Picard-Lyapunov type, we next integrate (12.61) three times with conditions (12.81) to obtain 
the equivalent integral equation 

(2-9) f(y) = 1 - C(y - yf - J J j y y ° h(f(y),y) (dyf. 

It is convenient to rewrite (I2.9P for the functions v(y) given by shifting 
(2.10) f(y) = l-C(y -y) 2 + v(y). 

Then the integral equation reads 

(2.H) v(y) = M(v)(y) = - J J J y yo h(v(y),y) (dyf, 

where we denote h(v, y) = h(l — C(y — y) 2 + v,y). 

We next introduce a functional framework, which is suitable for integral equation with 
singularities; see (33J, § 5.8], where these techniques were applied to uniform stability 
problems for degenerate Hamilton- Jacobi and singular perturbed parabolic PDEs. Thus, 
let us fix a small 5 > and set Ig = [y — 5, y ). We consider the integral equation (12.111) 
in the functional space 

(2.12) X p = {u6 C{I s ), v(y ) = v'(y ) = 0, pv G L°°(I S )}, where p = 

We endow X p with the following sup-norm: 

(2-13) \v\p = SMV yeh {p(y)\v(y)\}, 

i.e., we introduce the singular weight fa~^3 ~~ > 00 as y ~ *■ ^ ne standard framework 

of the space C(I$) of continuous functions with necessary boundary conditions at y = 0. 
With this distance (12.131) . X p becomes a complete metric space being a closed subspace 
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of the Banach space C P (I$) of twice differentiable functions at y — yo, for which the 
semi-norm 

(2.14) \v\ p = sup yeIs {p(y) \v(y) - v(y ) - v'(y )(y - y ) - \ v"(y )(y - y f\ } 
is finite. The space C p (Is) has the natural norm 

Hp = HVo)\ + \ v '(yo)\ + W'(Vo)\ + \v\ P . 

Firstly, it is an easy exercise to see that two terms in (12.61) being rewritten for v, 
immediately yield 

(2.15) M : X p -> X p for all small 5 > 0. 

Indeed, (12.61) implies that, for y « y , h ~ ^ + O((y — y)), so, on triple integration in 
(12.111) . we obtain the result: 

0<M(v)< [f + o(l)](y -yf<y (y -yf in I s . 

Secondly, checking the contractivity of A4 on X p , one can see that the first term h± in 
(12. 6p does this as inducing non-singular and analytic operator. Consider the second term 
h 2 therein. Taking arbitrary v\ j2 G X p and denoting by fi t2 their /-counterparts due to 
the change in (I2.10p . setting for convenience Av — v± — v 2 — fi — f 2 will provide us by 
standard manipulations similar to the Lagrange formula of finite increments the following: 

M/i) ~ M/2) = h 2 {vi) - h 2 (v 2 ) 
_ y , J(i-/i) _ y _ /(1-/2) 

4(1+/!) "T 4(1-/?) 4(l+/ 2 ) 4(1-/|) 

t 2 ' 16 ) y(/2-/i) 1 K1-/2) 1 /(1-/2) J(i-/ 3 ) 

(l+/l)(l+/ 2 ) ^ 4(1-/?) 4(1-/?) f 4(1-/?) 4(1-/2) 
_ 2/Ai; /Ap J(l-/ 2 ) ( 1 1 \ 

(l+/i)(l+/ 2 ) W^ff)^ 4 li=/f T=7fJ- 
Obviously, for small 5 > 0, we may assume that, in Is, 

l<l + /i, 2 <2, I^^^l, 



(2.17) c 



(yo - yf < 1 - h,2 < 2C(y - yf, < 



2 vyu _= - JV vyu , C ( yo -y) 2 - (l-/» a ) - C( yo - y y 

Then, denoting here and later on by A > various constants that are independent of 5, 
we obtain the following estimate: 

- kfoOl < f + / |A,| + I Sf* 

<A\\Av\+ r -^J\Av\ + ^]. 

— LI I (yo-y) J 1 1 2/0-2/ -I 

Using the fact that, in the metric (12.131) . 

(2.19) \Av(y)\ < \Av\ p (y Q -yf, 

we obtain that 



(2.20) \h 2 ( Vl ) - h 2 (v 2 )\ < A\Av\ p [(y - yf + (y - yf] < A(y - yf\Av 
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Finally, using the metric (12.131) again and substituting (I2.20p into (12.111) . we then obtain 
the following principal bound: 

\M( Vl ) - M(v 2 )\ p < p(y) J J /;° A(y - yf |A«|„ {dyf 

(2 ' 21) < \ Av \p a j^w fff?iv°-v) 2 W <£s 2 \Av\ p . 

This implies that 

(2.22) M is a contraction in X p for all small S > such that ^ S 2 < 1 in (12.211) . 

Hence, by Banach's Contraction Principle (see e.g. [351 p. 206]), for all sufficiently small 
5 > 0, (12.151) and (12.221) guarantee existence and uniqueness of a suitable local solution of 
(I2~TTD in X p for any fixed C > 0. □ 

2.2. Analyticity at y = y$. This is also a principal question of TFE theory. We begin 
with two simpler illustrations. 

1. Explicit analyticity for the standard TFE. Consider the "standard" TFE 

(2.23) u t = -{uu xxx ) x in KxR + , 

which has the same type of degeneracy as (11.131) but now at the equilibrium {u = 0}, so 
that equations near degeneracy sets are approximately connected by the change u \— > 1—u. 
Looking for the source-type solution of f!2.23j) leads to a simpler ODE in divergence form: 

(2.24) u s (x,t)=t-*f(y), y=-^ => -(//"')' + I = 0. 

The ODE in (12.241) near the interface is again approximately connected with (12.11) by 
f i — > 1 — f. Integrating (12. 24f) ones and dividing by / > yields the following explicit 
solution detected by Smyth and Hill in 1988 [B]: 

(2.25) -/r + |/y = =► f'" = \y, i.e., f(y) = JL (yg - y*)\ 

Thus, we see a finite polynomial solution (I2.25p . which is a perfect and simplest analytic 
solution at the interface y = yo, which, surely, is a singularity degeneracy point for the 
ODE in (12.241) . As an easy illustration, we observe that the interface at y = yo occurs via 
the geometric configuration when 

(2.26) an analytic function f(y) "touches down" the level {/ = 0} at some y = y > 0. 

Being formally extended by f(y) = for \y\ > y$ (in fact, this makes no particular 
sense for the FBP posed in {\y\ < y } only), the solution (12.241) solves another singular 
"Riemann problem" with a measure as initial data, 

u s (x, t) — >• c d(x) as t — >• + (c = ^j-) 

in the sense of distributions. Hence, (12.241) is a formal fundamental solution of the FBP. 

Unfortunately, as we have demonstrated, our ODE (12. ip is more complicated and does 
not admit an explicit integration. However, we are going to show that, for the FBP, 
analyticity of the similarity profile f(y) as in (I2.25P is not exceptional and exclusive, so 
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f{y) of (12.11) remains analytic and actually is close to that in (I2.25P as y — > y , i.e., with 
the change / i— > 1 — /, the FBP profiles for (12.11) for y w y 

(2.27) are given by small analytic perturbations of the analytic solution in (12.251) . 

Therefore, the explicit finite polynomial FBP solution in (12.251) is a good illustration of 
the origin of the real analyticity of other non-explicitly given solutions. In its turn, the 
analyticity has also a strongly motivated natural origin in complex variable presentation: 
2. Analyticity for complex-valued ODEs and Weierstrass theorem. The origin of the an- 
alyticity of the FBP profiles can be also seen from the ODE theory in the complex plane 
z = x + \y G C, where x now stands for y in (12. ip . Namely, consider an analytic function 

(2.28) / = f(z) = u(x, y) + iv(x, y) (u x = v y , u y = -v x ) 
satisfying the complex ODE 

(2.29) -(^-f 2 )fz ZZ )z + \zf z = in C. 

Since (12.291) contains analytic nonlinearities only, it admits a family of analytic solutions 
in a neighbourhood O of the origin z = 

(2-30) f(z) = E (fe > ) tkz\ 

with the following relation on the expansion coefficients: 

c k+4 (k + l)(fc + 2)(fc + 3)(fc + 4) 

(2 ' 31) - E c i c k+ ^i- j c j l(l-l)(l-2)(k + l) = lkc k , k>0. 

l>3,0<j<k+4-l 

Choosing an analytic solution that is better associated with the FBP expansion (I2.2p . i.e., 
Co = 1 and C\ = 0, 

(2.32) f(z) = 1 + E(fc> 2 ) °k( z ~ z o) k , where z = (x , 0), 

yields a similar algebraic system, 

E c j c l c k+4 + j l(l-l)(l-2)(k + l) 
(2 33) '>3, i<j<fc+4-/ 

= | kc k + \ x (k + l)c fc+ i, k > I. 

In fact, the analyticity of the solutions, i.e., non-zero radius of convergence of the series 
in (12.301) and (12.321) can be seen from the infinite algebraic systems (I2.3ip and (I2.33p . In 
other words, one needs to show that the recurrent relations generated by (I2.3ip or (12.331) 
admit solutions {c k } without huge growth such as 

(2.34) limsup^^ \c k \k = oo (then R = 0). 

A direct proving that (I2.34p does not take place for the above power series is not straight- 
forward (but indeed seems doable), so we will follow a simpler approach. 

Meantime, we note that, on the other hand, iterating the integral equation equivalent 
to (I2.29P via the simple iteration starting with an analytic initial data fo{z) in O (this 
will be done shortly for the equivalent integral equation in the real case, so we do not 
discuss details) will give us a sequence of analytic functions {fk(z)}. In case of the uniform 
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convergence of the series along a subsequence in a neighbourhood of z = (again, the 
convergence is naturally proved for the integral counterpart of the equation), 

fk{z) -> f{z) as k = kj -> oo, 

will provide us with an analytic solution, as a manifestation of Weierstrass' classic theorem 
(the limit of a uniformly converging sequence of analytic functions in O is also analytic). 

This is the origin of analytic solutions of the degenerate ODE (12.291) in the complex 
plane. Let us trace out the link to the real case. In the variables in (I2.28p . currently using 
the differentiation f z = u x + iv x , the ODE (12.291) reduces to a system for Re / and Im/, 

J-[(l - u 2 + v 2 )u"' + 2uvv"')' + \ (xu' - yv') = 0, 

( "" >J) {-[(1 - u 2 + v 2 )v"' - 2uvu"'}' + \ {yv' + xv') = 0, 

where ' = D x . It then follows that: 

(2.36) (J2HD is (E29D at the axis {y = 0}, on which v{x, 0) = 0. 

Indeed, then also v x (x, 0) = u xx (a;, 0) = ... = that make the second equation in (12.351) 
tautological, and then the first one for f(x) = u(x, 0) coincides with (12.11) . where y i— > x. 
One can see that an analytic solution of (12. 29f) satisfying (12.361) must have the form (12.301) 
with all real expansion coefficients^ {c^}. Then, among those analytic solutions, there 
can be some, which we have introduced in (I2.32p . corresponding to a "touching down 
configuration" of u(x,0) as in (I2.26p . A proper solvability of the analytic elliptic system 
(12.351) with the condition at the x-axis in (I2.36j) . i.e., in the class of R-analytic solutions, 
and with the geometry as in (12 . 26[) is not a part of the present study. However, this shows 
a key analyticity link between complex and real ODEs and the role of Weierstrass' theorem 
in constructing analytic solutions of such degenerate ODEs with analytic nonlinearities. 
3. Main result . Thus, we now turn our attention to the real valued ODE (12.11) and prove: 

Proposition 2.2. The local solution of Proposition [Ol is real analytic at y = y . 

Proof. Uniqueness. First, let us note again that in (12.61) . the main first non-singular 
term itself produces existence of an analytic solution. At the same time, the second 
perturbation term is also analytic in the sense that it maps admissible (i.e., those that 
do not create singularities in the integral) analytic functions into analytic. Therefore, the 
ODE (12. 6p admits a unique formal analytic expansion. Indeed, setting 

(2.37) /(*/) = 1 + E (i > 2 ) 9 A ^ wh ere A = y - y (c 2 = -C < 0) 
and substituting this formal expansion into (12. 6p yields another expansion 
(2-38) h(f(y), y) = 4(w } (1+/) ft f = E (fc >i) b k A k , 

with a uniquely solvable relation between coefficients {c^} and {6^} as in (I2.3ip . Inte- 
grating (12.61) three times leading to (12.111) preserves this formal power expansion. Thus, 
this yields the uniqueness: if an analytic solution (I2.37P of (12.111) exists, it is unique. 

3 It is convenient to call such solutions R-analytic, with "R" standing for "Real coefficients" (not to 
confuse with "real analytic" to be treated below). 
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Note also that this implies that the solution constructed in Proposition 12.11 is at least 
C°°, i.e., applying finite expansion analysis with the remainders of order 0(A l ), with 
arbitrarily large / 3> 1, yields that there exist all the derivatives 

(2.39) 3f^(y )=j\ Cj for all j > 0. 

Of course, this is easily seen from (12. 6p by differentiating as many times as necessary: 

f"(y ) = -2C, f"(y ) = y f, / (4 %o)=0, = -gg, ... , 

which lead to the expansion (12. 2ft . 

Existence . Let us show convergence of the series (12.371) . First of all, there exists an a 
priori bound for (|2.6|) : for any suitable solutions / such that v 6 X p , in 7^, 

(2.40) < h(f, y) < a + / 2C(y - I/) 2 dy < f + «p < y . 
We now perform an analytic simple iteration of the integral equation (12.111) . 

(2.41) v k+1 = M{v k ), k = 0,1,2,..., 
by taking analytic initial data 

(2.42) Vo = OeX p , i.e., f (y) = 1 - C(y - y) 2 . 

By Banach's Contraction Principle, we then obtain a sequence of functions {fk(y)} such 
that the corresponding polynomial sequence {vk(y)} given by (12. 101) satisfies as k — > oo 

(2.43) ^fc(y) — * v(y) in X p and uniformly in Is. 
These functions are given by converging power series 

(2-44) /fe(2/) = l + E > 2 )4 Ai - 

By ( I2.40p . the polynomial sequence {fk(y)} is uniformly bounded in 1$. 

Finally, we then need to prove that the radiuses of convergence of the approximating 
series ( I2.44p are bounded from below: 

(2.45) R k > 5 > for all k > 1. 

To this end we need some more calculus concerning the integral operator involved. 
Thus, according to (I2.9p . by the change 

( = y o -y>0, f(y) = l-C( 2 (l + e(()) 

equation (12.61) reduces to 

(2.46) C(C 2 e)'" = h(e, () = ^ + 4(1+ / )(1+£) ^ J ? C 2 (l + e)d(, 

where now the right-hand sided is analytic in e but is not regular in ( in the last term 
with formal singularity ^ , which we have to pay the main and the only our attention. For 
simplicity and without fear of confusion, we still keep / in the terms that are obviously 
regular analytic and provide no difficulties. Then changing the independent variable in 
the left-hand side 

(2.47) C = fa) > =^ (C 2 e) w = CV" + ... = e»> + ... 
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we set <$' = 2 / 3 , i.e., Q = 4>{s) = (|) 3 gives an analytic change of variables, to obtain 
(2-48) Ce»> + ... = h(e, (f ) 3 ) = ... + i(TT7 ^ Tiy £ £ , 8 (1 + e) ds, 

where we omit lower-order regular terms in s, which supply us with guaranteed analytic 
s-expansions. Actually, the integral representation (I2.48P actually completes the proof of 
analyticity. Indeed, we already know that the solutions is unique, and then substitution 
into (I2.48P the analytic series gives a rather standard relation between expansion coeffi- 
cients that guarantee the non-zero radius of convergence. Let us present some comments. 

Consider the iteration (I2.4ip in the e-variable that gives us the required sequence 
{^k, k > 0}, which thus is given by the recursion 

(2.49) Ce»> +1 + ... = h(e k , (§) 3 ) = ... + 4(1+/ ^ (1+£fc) £ £ * 8 (1 + e k ) ds. 

Let us perform the first step to estimate the radius R\ > 0. Let C = yo = 1 and take 

£o = and /o(C) = 1 - C 2 = 1 - (f) 6 - 
This yields that £i(s) is given by the triple integration of the equation 

f'" + - I 1 1 f s n 8 dn 

fc 1 t- ... — ... t- 2 -(|)6 s 6 Jo 6 Ui ' 

where we again keep the "non- analytic" term only that, by assumption, may eventually 
vanish the radius of convergence of approximations £fc(s). We also omit all the constants 
that do not affect this radius. Thus, E\ is given by integrating 

(2.50) + ... = ... + 2r§)6 #P = 3 • 2« > 1, 

where we indicate the radius of convergence for the power series for £i(s), which is deter- 
mined by the formally non-analytic term (the actual radius can be smaller if the analytic 
terms define a smaller radios, that does not happen). Thus, £i(s) admits the unique 
analytic continuation denoted again by Si(z) into the disc B\ = {\z\ < Ri} C C. 

Consider the transition k h- > k + 1 via (12.491) . Let s^(z) be analytic in 5^ = {|z| < Rk}, 
where Rk > 0. By the definition of e, we may assume by (12.431) that 

(2.51) £fc(0) = and |ejfc(s)| < \ in 7<j (actually |£fc(s)| < \s\), 

and moreover we also may assume that the analytic continuation £k{z) also satisfies this 
estimate in Bk- Consider the right-hand side of ( 12.49H . Obviously both functions 1 + fk 
and 1 + Ek are analytic in B^, and so do y^r- and j^- provided that these do not vanish 
in Bk that is true by (I2.5ip for small 5 > 0. The integration in (I2.49p also does not change 
the radius of convergence. Then, overall, the product of three analytic functions in (12.491) 
preserves the same radius of convergence Rk, so that by triple integration £ fc+1 (s) will 
inherit at least the same uniform estimate of i2£7* from below. 

Thus, inside this interval of uniform convergence of all the series, we get the desired 
result: as k — > oo, uniformly in Is, 

(2.52) fk(y) = 1 + E (i > 2) S* Ai - 1 + Ey> a) ^ = f{v), 
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where f(y) is the unique solution from Proposition 12. II This implies that f(y) is analytic 
at the degeneracy point y = y . □ 

It is worth mentioning that in the above analyticity proof, we essentially used the 
existence-uniqueness via contractions in C p established in Proposition 12.11 On the other 
hand, the presented estimates of the radiuses of convergence of {ek(z)} can also provide 
existence and uniqueness of analytic solutions, but this issue should be pre-converted into 
the complex variable nature for the ODE (12.291) . Then, this analysis will precisely em- 
phasize the application of the already announced Weierstrass theorem. However, proving 
non-existence of non-analytic solution profiles will anyway require a contractivity study 
of the integral equation such as (12. Ill) in metric spaces of functions of finite regularity. 

2.3. First similarity profile: numerics. Concerning existence of finite interface in 
the degenerate ODE such as (11.171) . we refer to first results in [6j [8]. Existence (and 
uniqueness) of the first similarity profile / = fi(y) G [0, 1] depends on a 2D-2D matching 
argument, where, using two parameters in the bundle (12.31) . we need to satisfy three 
conditions (I1.19P with the unknown interface point yo (a free parameter); see below. 

Shown in Figure [1] are results of numerical shooting for the FBP (I1.17I) - (I1.19I) . We 
indicate here a number of profiles (dotted lines) satisfying the anti-symmetry conditions 
at the origin (11.181) and two conditions at the right-end point 

(2.53) /(x/o) = 1 and f"'(y ) = 0. 
We next vary the length y > to get the zero contact angle 

(2.54) f{y ) = 0, 

and then, by construction, all three conditions (I1.19P are valid. We obtain the first 
profile fi{y) (the boldface line) satisfying the necessary range condition (11.231) . which was 
numerically clearly unique, with the interface at 

(2.55) j/oi ~ 2.35. 

The above shooting procedure also indicates that there exists the unique second FBP 
profile f2(y) with the interface at 

2/02 ~ 3.75. 

We return to this multiplicity problem of the FBP profiles in Section [3] after introducing 
and studying the unique oscillatory similarity profiles of the Cauchy problem, to which 
FBP ones turn out to have a direct relation. 

2.4. On existence and uniqueness of the first monotone FBP j\. We begin by 
mentioning that the mathematical techniques developed in Bernis [6\ and Bernis-McLeod 
[S] as well as in [TUl EE] , are essentially directed to third-order ODEs obtained on integra- 
tion ones via conservation laws, so these do not directly apply to the truly fourth-order 
ODE (I1.17p . Nevertheless, the methods and results therein can be useful. In particular, 
Bernis energy estimates [H § 7] settles finite propagation property for such ODEs. 

Shown in Figure [1] is the actual mathematical strategy to obtain the necessary profile 
fi(y) to get the following: 
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Shooting FBP profiles for n=1, m=2 



f(y)i 




FIGURE 1. Construction of the unique first similarity profile /i(y) for the FBP. 

Proposition 2.3. The problem ( 1.17 )- (\1.19t\ admits a solution fi(y) satisfying 

(2.56) fi(y)>0,fi<0 on (0,y ), so < h{y) < 1. 

Proof. Thus, we shoot from the unknown interface point y = y$ > according to the 
bundle ( 12. 2p . and denote the solution as 



(2.57) f = f(y;yo,c), 

which can be continued until the regular point y = 0. One can see that the ODE (11.171) 
does not admit local singularities at finite points. Our further analysis uses the continuous 
dependence of f(y) on the parameters yo and C. Note that, as we have shown, the 
dependence is also continuous in the expansion ( 12. 2ft near singular (analytic) points. 

We next need a few simple observations concerning such shooting orbits: 

(i) We want to know that an inflection point of f(y), where 

(2.58) f"(yi)=o, yi e[o,y ), /(j/i)e[0,l), 

actually occurs. This is seen from the ODE and, moreover, from the expansion (12.21) that 
is obtained from it. Indeed, we have that, for y rs y$ , 

(2.59) f"(y) = -2C+ly {y -y) + ... = at y x = y (l - ^ + ...), c 
The second condition in (I2.58P then demands 



^< 1. 

Vo 
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Existence of fi(y) is proved if, for some (y ,C), 

(2.60) Vl = 0. 

(ii) We next use again the easy fact that follows from the ODE ( 11.171) that consists of 
two terms only with clear depending sign properties. Thus, ( 12. 4ft implies the positivity, 

(2.61) -((i-f 2 )n(y) = Uy V0 f'y>^ 

in a neighbourhood of y^, where f(y) is strictly increasing. Therefore, /'" < there, so 
that f"(y) is strictly decreasing while f'(y) > 0. Hence, we conclude that a monotone 
function f(y) cannot have more than a uniquJI inflection point y\ = yo(yo,C), at which 
(12.581) holds, in the connected domain of monotonicity of f(y). 

(iii) It is easy to see from (12.21) that, for any fixed y > 0, 

(2.62) /" < on [0,y Q ) provided that C > y*, 

so that the inflection point in (I2.58P disappears for sufficiently large C. It is also easy to 
see that the type of disappearance of the point I(y ,C) = (yi(yo, C), f(yi(yo, C))) from 
the domain Si = {y > 0, / G (0, 1)} is different for y ^> 1 and yo -C 1. Namely, for 
yo 3> 1, the point I(yo) passes through the boundary axis {/ = 0} as C increases, while 
for y «C 1, this is done through the boundary segment on {y = 0}. 

Using the above estimates, by continuous dependence on parameters y and C of this 
analytic family {/}, we conclude that there exists a yo > such that (I2.60P holds and the 
fi(y) satisfies (I23BD . □ 

We expect that, in similar lines, there exists a proof of existence and uniqueness for 
higher-order TFEs leading to ODEs such as (11.321) that still contain two operators with 
related "signs" . Note that the total number of parameters increases with the order 2m so 
the shooting approach becomes more and more involved. For instance, for m = 3, i.e., for 
the sixth-order TFE (11.301) . we then have a three-parameter shooting, with parameters 
/'(0), /'"(0), and/( 5 )(0). Nevertheless, the uniqueness of fi is surely observed in numerical 
experiments in Section [HI 

2.5. On stability of the similarity solution in PDE setting. We now show how to 
check whether the similarity solution (I1.16P with / = fi(y) is stable for the TFE (11.131) 
with respect to small perturbations of the data. We introduce the rescaled variables 

(2.63) u(x,t)=v(y,r), y = fji , r = lnt, 
where v solves the equation with the operator B from (I1.17p . 

(2.64) v T = B{v). 

Since r — > — oo as t — > 0, where u(x, t) approaches a slightly perturbed data S + (x), we 
need to show that these perturbations are exponentially decaying as r increases. Thus, 



It seems that this is key for the uniqueness, but we still cannot complete the argument, where, as 
usual, an extra monotonicity-\ike property for the ODE is necessary. 
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one needs to check the spectrum of the linearized operator 

(2.65) B'(f)Y = -((1 - f 2 )Y"J + (| + 2ff")Y' + 2(f f")% 

that is posed on y G (0, yo) with the anti-symmetry regular conditions at y = 0, 

(2.66) F(0) = y"(0) = 
and the free-boundary conditions 

(2.67) Y = Y' = at y = y . 

For simplicity of linear stability analysis, we assume that the free boundary is fixed at 
= yo, so we exclude variation of the free boundary that are not principal. 

Operator (I2.65j) is not symmetric in any space such as L 2 , so it does not admit a self- 
adjoint extension. It is rather difficult degenerate singular operator that at the singular 
point y = yo has the principle part with a quadratic degeneracy according to (I2.2p . 

B>(f)Y=-2C((yo-y) 2 Y">y + .... 

Therefore, solving the principle part of the eigenfunction equation, 

(2.68) - 2C((y - yfY'")' + ... = XY « 0, 

one obtains that the strongest singularity admitted by (12.681) is bounded, 

(2.69) Y s (y) ~ (y -y)ln(y -y) + ... for y^y - 

Obviously, the singular part (12.691) does not satisfy (12.671) . 

Hence, the conditions (12.671) themselves define an operator extension with a discrete 
spectrum. This follows from classic theory of ordinary differential operators [36J, and 
also can be easily seen geometrically. Indeed, looking for eigenfunctions solving (12.681) 
is equivalent to shooting from the point y = 0, where in view of (I2.67P we are left with 
two parameters Y'(0) and Y"'(0), and actually, since by linear scaling we can always 
put Y'(0) = 1 (if it is not zero, which is a special non-generic case). Therefore, the 
parameter Y"'(0) G C together with A G C are, in general, sufficient to satisfy precisely 
two conditions (I2.67P at y = yo- This gives an analytic system for A |36j, and shows that 
a compact resolvent in a suitable functional space as a meromorphic function is detected. 

We restrict our attention to a necessary bound of the spectrum of (I2.65p . 

(2.70) B'(/)y = XY. 

As usual, multiplying this by the complex conjugate Y in L 2 , taking the conjugate and 
multiplying by Y, and summing up both yields 

(2.71) - f(l - f 2 ) \Y'f + fl-(ff')'} \Yf + J [(ff"J ~ |] \Y\ 2 = / \Y\ 2 . 

It is crucial that the first, highest-order term is negative, so the rest of the positive terms 
are expected to be estimated via interpolation. Actually, there are only two positive 
components: (i) in the second term: 

(2.72) - (//')' = -//' - (f)\ where - //" > by ([236}, 
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and (ii) in the last term 

(2.73) (//"')' = // (4) + /'/'", where f f" < but the sign of ff^ is unknown. 

Note that close to y = yo, the interpolation is possible regardless the degeneracy of the 
principle part in (12.711) . 

Thus we conclude from (12.711) that 

(2.74) Re A < 0, 

provided that a suitable (Friedrichs) self-adjoint extension of the following symmetric 
operator satisfies (a proper setting is standard, [36]): 

(2.75) Q P (f)Y = -[(1 - f)Y"\" + [(ff)'Y']' + [(ff"J - J]K < 0. 

For the first FBP profile fi(y), which is not given explicitly, (12.751) can be checked numer- 
ically and is rather plausible. At least, even without careful and not that easy spectral 
numerics, we observe that the self-adjoint operator in (12.751) exhibits a clear tendency to 
be negatively determined for reasonable FBP profiles / = fi(y) from (12.561) . 

3. The Cauchy setting for RP-1: oscillations and similarity profiles 

We begin by noting that the CP setting for the TFE (11.121) assumes no analyticity 
unlike the TFE setting discussed in Section [21 For the CP, we observe oscillatory and 
changing sign behaviour near the interfaces. 

3.1. Local solutions of changing sign. The oscillatory local behaviour for (I1.17P is 

more complicated. Here we follow [24]. We introduce the oscillatory component <p(s) as: 

(3.1) f(y) = l-(y - s = \n(y - y), 
where ip(s) solves an autonomous ODE with a discontinuous nonlinearity, 

(3.2) P 3 ((f) = ip'" + 6</ + IV + Qip = -A sign;/?, where A = \ y > 0. 

We look for a periodic solution of (13.21) . which by (13.10 describes the oscillatory nature of 
solutions in the CP. Let us first note that technique based on calculation of rotation of 
vector fields for the dynamical system (13.21) does not apply here; see [351 P- 45-53]. In this 
case, existence of multiple periodic solutions depends on the careful analysis of the index 
of a guiding function V (it always exists and can be calculated explicitly by using the 
linear form of the operator at infinity). Indices of such periodic solutions are unknown, 
the main results in [221 p. 52] seem not applicable for present problem. 

Therefore, a direct shooting approach is effective here. We consider a 2D bundle of 
orbits <f(s, Ci, Cz) satisfying the conditions 

(3.3) ?(0) = 0, </(P) = C 1 , </(0) = C 2 , 

where G K are arbitrary parameters. Let us state some key properties implying 
existence of a periodic orbit. 
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(i) (p(s, Ci,C 2 ) is uniformly bounded for all s > 0. Obviously, in view of boundedness 
of the nonlinearity, (p(s, Ci, C 2 ) is locally well denned. If, on the contrary, (p(s, C\, C 2 ) is 
unbounded, we should have that the linear counterpart, 

</' + 6y/' + lV + 6y? = 0, 

must admit unbounded solutions. Setting ip(s) = e MS gives the characteristic equation 

$0) = fi 3 + 6/i 2 + + 6 = 0. 

Since, $'(//) > for /i > 0, $(0) = 6, and 

$'(>)= for /i = fi± = |(6 ± v 7 ^), with $(// ± ) = ^0.3849... , 

we have that all three roots of are negative. So the solutions cannot grow without 
bound as s ^ +oo. 

Actually, this means that, as a consequence, that the fourth-order ODE (13.21) with 
bounded coefficients is a dissipative DS having a bounded absorbing set. Dissipative DSs 
are known to admit periodic solutions in a rather general setting [331 § 39] provided these 
are non-autonomous (so the period is fixed). For the autonomous system (13.21) . the proof 
in [2U § 7.1] can be completed by shooting. 

(ii) ip(s,Ci,C2) is oscillatory as s — > +oo, i.e., has infinite number of zeros in any 
neighbourhood of s = +oo. Indeed, if, say, (p(s, C\, C2) > for all s ^> 1, we will have 
there the linear ODE 

(ff" + Q<p" + IV + Qif = -1, 
admitting the particular equilibrium solution <p = — ~. Hence, by the above stability of the 
orbits of the homogeneous linear equation in (i), we conclude that <p(s, G\, C 2 ) ~ — | < 
for all s ^> 1, from whence comes a contradiction. 
We now state the main result. 

Proposition 3.1. (i) Equation (\3.2l) admits a unique non-trivial periodic solution <fo(s) 
of changing sign, which is stable as s — > +00; and 

(ii) The periodic orbit (po(s) of (\3.2$ is the unique bounded connection with s = —00. 

Proof, (i) Existence and uniqueness have a pure algebraic proof, [2H § 7.4]. Stability 
follows from the above analysis of the linearized problem, (ii) This also implies that the 
stable manifold of (fo(s) as s — » —00 is empty, so, by the geometric analysis in [23], <^o( s ) 
is the only bounded orbit connecting s = —00. □ 

Figure [2] shows that the periodic solution (po(s) is stable, up to translations in s, as 
usual, for the ODE (13.21) as s — » +00. 

In view of (ii), the expansion (13. ip with (p = ifo(s + Sq), Sq G K. is arbitrary, is the only 
maximal regularity connection with f(y) = for y > y . Therefore, this describes the 
generic robust structure of the multiple zero at the interface of arbitrary similarity solu- 
tions of the TFE-4 under consideration. We expect that this structure remains unchanged 
for general solutions of the PDE (I1.12p . which is a difficult open problem. Observe that 
according to (13.11) . the actual regularity of f(y) is C 2,1 that is better than for the above 
FBP with the less smooth C ' l behaviour in (12.21) (if we formally set f(y) = l for y > yo). 
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FIGURE 2. On the stable periodic orbit of (pH|) . with y = 4. 

As we have mentioned, the regularity in 03. ip can be attributed to the Cauchy problems 
and in fact is the maximal regularity which is admitted by the ODE under consideration, 

3.2. Similarity profile fo(y) for the CP. We begin by noting that finite propagation 
in the ODEs such as (I1.17I) . i.e., existence of a finite interface yo, is well-known for a long 
time; see first ODE proofs in [HIE], and more general energy estimates for PDEs in [51 |4"T] 
and survey in |31j . 

As in the FBP, the bundle of asymptotics (13. ip is 2D comprising parameter {y , s }, 
where so is the translational constant in the oscillatory component, <fo(s) i— > (fo(s + sq)- 
Of course, the matching problem for the CP with the oscillatory bundle (13. ip gets more 
difficult. 

In Figure [3], we show the unique similarity profile fo(y) corresponding to the CP. Figure 
H] explains the character of oscillations of the f(y) about / = 1. This behaviour well- 
corresponds to the expansion (13.11) . This shows a rough estimate of the interface location 

(3.4) y ,CP ~ 6.7, 

which is difficult to improve numerically in view of the oscillatory behaviour, the necessity 
of regularization parameters, and sufficiently low tolerance of convergence that was applied 
to the dynamical system (equivalent to (11.171) ) with the essentially non-symmetric matrix. 

3.3. On a countable set of FBP profiles. Figured] indicates an interesting "geometric" 
relation between the first FBP profile and the CP one (n = 1). For further convenience, 
in Figure \5[ we present the enlarged version of this figure with extra details, where we 
include shooting of the third FBP profile fz{y) with interface positions at, respectively, 

y i ~ 2.35, T/02 ~ 3.75, and y 03 ~ 4.74, 

which all are smaller than the CP interface location (13.41) . By boldface dotted line therein, 
we denote the CP profile fo{y)- 
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The Cauchy problem 




FIGURE 3. The CP profile f (y) satisfying (fTTTD . (TOTl) . 



The CP: oscillations about f(y) = 1 



The CP: oscillations about f(y) = 1 , enlarged 





(a) oscillations, enlarged 



(b) oscillations, further enlarged 



FIGURE 4. Enlarged oscillations of /o(y) from Figure [3] about the equilibrium 
/ = !■ 

Namely, we see that the touching point of fi{y) at the level 1 belongs to the interval of 
the first oscillation of the profile for the CP, and the same is true for f 2 and /3 in Figure 
Moreover, we expect a further extension of this geometric property (an open problem). 

Conjecture [3J.2. The FBP ( | 1. 17 )- (\ 1. 1% admits a countable family of different solutions 
{fkiy)}, where fiiy) G [—1, 1], which converges to the CP profile: 

(3.5) fkiy) -> f (y) as k -> oo 

uniformly on any interval from [0,?/o,cp)- 

In Figure [6], for clarity, we present a formal scheme of this geometric interaction between 
FBP and CP profiles, which partially has been reflected in Figure [51 This was already 
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FBP and CP profiles for n=1 , m=2 




0.995 



FIGURE 5. The CP f {y) and first three FBP profiles for n = 1. 

seen in Figures [TJ El and will be clearer explained in Figure [14] obtained numerically for 
the TFE-6 (m = 3 in (11.301) ) admitting more oscillatory patterns. In particular, these 
figures show that the interface points y k of the FBP profiles satisfy 

(3.6) y ok -> y CP as k -> oo. 

Recall that any FBP profile fk{y) being put into the similarity formula (11.161) gives a 
solution of the problem (I1.12p . (11.141) . Therefore, in the FBP setting, Riemann's problem 
for the TFE-Jj. has an infinite countable number of solutions (and only the first profile 
fi{y) has the physical range [—1,1]), while for the Cauchy setting, such a solution is 
expected to be unique (and violates the range [—1, 1]). 

4. On extension at the branching point n = in the Cauchy problem 

Here, we demonstrate another approach to Riemann's problem in the case of the Cauchy 
problem for the TFE-4 (11.221) with parameter n > 0. Then the similarity solution has the 
same form (I1.16p . there f(y) solves the ODE 

B n (/) = -(|l-m"') , + \fy = for y>0, 

/(0) = /"(0) = 0, /(+oo) = l. 

By our local analysis, we expect the solution with the finite interface to have a sufficiently 
regular connection with the equilibrium / = 1, as explained in Section [31 

Our basic idea is to show that the solution of (14. ip can be extended from the solution 
/o in (11.281) of the linear problem for n = 0, which then becomes a branching point for 
(14. ip ; see [351 P- 371]. To justify this, we will need some extra spectral properties of the 
linear operators involved. 
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FIGURE 6. An extra schematic explanation of a countable set {fk(y)} of FBP 
profiles approaching the CP one (the bold line); cf. Figure 

4.1. Similarity profile and spectral properties of B for n = 0. Thus, as we have 
seen, for n = 0, the problem ( 14. II) in the CP setting takes the form 



where B is the linear operator in (11.261) that defines the rescaled kernel of the fundamental 
solution of the parabolic operator D t + D\. 

The necessary spectral properties of the linear non self-adjoint operator B introduced 
in (ll.26p and the corresponding adjoint operator B* are of importance in the asymptotic 
analysis, as explained in [20] for similar general 2mth-order operators (see also [231 § 4]). 
In particular, B is naturally defined in a weighted space L 2 p (W), with 



The corresponding eigenfunctions are normalized derivatives of the rescaled kernel, 



B (f) = -f^ + {f'y = Q for y > 0, 
/(0) = /"(0) = 0, /(+oo) = 1, 




p{y) = e a ' y ' , where a £ (0, 2d) is a constant, 




(4.5) 



My) 




FW(y), 1 = 0,1,2,.... 



The adjoint operator 



(4.6) 



B* = -Dl - \ 
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has the same spectrum (14.41) and the polynomial eigenf unctions 
(4-7) #(v) = ^Efc? J JW 1 = 0, 1,2,..., 

which form a complete subset in L 2 pt (M.), where p* = -. Similarly, the domain of the 
bounded operator B* is the Sobolev space ifi(R). In particular, it is easy to compute 

(4.8) % = i, v? = v, r 2 = 73 y 2 > vs = 75 y 3 > ^ = 7= (</ 4 + 24), 

etc. As B, the adjoint operator B* has the compact resolvent (B* — A/) -1 in L^»(R). It is 
not difficult to see that, on integration by parts, the eigenfunctions H4.5|) are orthonormal 
to polynomial eigenfunctions {"0/*} of the adjoint operator B*, so 

(4.9) 

where (-, •) denotes the standard (dual) scalar product in L 2 (IR). 

Thus, in the necessary restriction to odd functions denoted by L 2 p {R + ), according to 
f)4.3p . we have that B has the discrete spectrum 

(4.10) <7(B„) = { - ^ < -|, / = 1,3,5,...} 

and a complete and closed set of eigenfunctions {ipi}- Hence, the linearized operator at 
n = is strictly negative in view of (14.101) that gives a good chance to extend the solution 
/o for n > small enough by applying classic branching theory in the study of the limit 
n -> 0+ in the problem (|Pjl : see 0SI P- 319] and pH p. 381]. 

4.2. Branching at n = 0. To this end, for n > 0, we write down (14.11) as a perturbation 
of the linear problem (14. 2p . 

(4.H) B f = g(f;n)= - P\ n - I) f"]' . 

Next, since B < is known to have compact resolvent in L 2 (M + ), we consider the 
equivalent integral equation for the function 

(4.12) Y = f-f , i.e., f = f + Y, 

which has the form 



n 



(4.13) Y = A(/ + Y; n) = B l g{f + Y 

Initially, we formally assume that, for small n < 1, the nonlinear operator on the right- 
hand side can be treated as a compact Hammerstein operator in L^-spaces, as classic 
theory suggests [35, § 17]. Recall that for n = 0, (14.131) is indeed a linear integral 
equation with a compact operator admitting the unique (up to a multiplier) solution / 
as in (OBI . 

Therefore, performing below necessary computations, we bear in mind a formal use 
of classic branching theory for compact integral operators as in (I4.13p . We also do not 
discuss here specific aspects of compact integral operators in weighted L p -spaces; see [20] 
for extra details. Note that, in a suitable metric, B is a sectorial operator, [27| § 5.1]. 
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Therefore, for derivation of branching equations, we use general branching theory in 
Banach spaces; see Vainberg-Trenogin [1HJ Ch. 7]. Beforehand, we need to discuss a few 
typical difficulties associated with the above eigenvalue problem. 

As (14.111) suggests, the crucial part is played by the nonlinearity 

(4.14) Q(f;n) = |1 - f 2 \ n - 1 in a neighbourhood of the point P = {f = l,n = + }. 

One can see that the derivative Q'f{f',n) is not continuous at Pq. Therefore, the standard 
implicit function (operator) theorem [l6j p. 319] does not apply. 

Nevertheless, the application of index-degree theory [35, p. 355] formally demands the 
differentiability only at the given point / = /o, n = 0, which is true (cf. more delicate 
computations below). Indeed, for n — 0, g(f,0) = 0, and the perturbation disappears in 
(14.131) at the branching point. 

On the other hand, in view of our difficulties with the regularity (and also with com- 
pactness of the operators involved), it is better to rely on Theorem 28.1 in [T3[ p. 381] 
that is formulated in the linearized setting, where the differentiability is "replaced" by 
the control of higher-order nonlinear terms as in (14. 14[) in a neighborhood. As usual, 
the key principle of branching is that the corresponding eigenvalue has odd multiplicity 
that can be easily checked in some cases (the even multiplicity case needs an additional 
treatment, which is also a routine procedure not to be treated here); see further comments 
below. Notice that the condition on the nonlinearity (14.141) is rather tricky to check that 
it satisfies the estimate (b) in Theorem 28.1 in [TJJ1 p. 381] for the variable x in / = F + x. 

Thus, we can proceed and conclude that (14.11) admits a solution f(y) for all sufficiently 
small n > and moreover these profiles form a continuous curve (an n-branch). 

We next detect a precise behaviour of this branch as n — ► + . 

4.3. Asymptotic expansion of the branch for small n > 0. We need to use in (14.111) 
the following expansion: 

(4.15) |1 - f\ n - 1 = ra In 1 1 - f 2 \ +o(n) as n -> 0+ 

which of course is not true uniformly on bounded intervals in / and should be understood 
in the weak sense. This is crucial for the equivalent integral equation (I4.13p . 

Proposition 4.1. For the function f — fo given by iH^M, in the sense of distributions 



and in the weak sense in L°°(WL + ) (and also in the sense of bounded measures in M + ) 

(4.16) i(|l-/ 2 | n -l) ^ln|l-/ 2 | as n^0 + . 

According to (I4.16p . analyzing the integral equation (14.21) . we can use the fact that, for 
any function <p G C, <p G L\R + ) (or (f> G C (M+)) 

(4.17) J (|1 - / 2 r - 1)0(2/) dy = n[f f(y) In |1 - f 2 (y)\ <j>(y) dy + o(l)] as n - 0+. 

Let us finish our formal computations that, for convenience, are performed for the 
differential equation (14. ip . Namely, substituting expansions (14.151) into (14. ip and using 
(14.121) yields the following perturbed equality: 

(4.18) B F = n(ln|l-/ 2 )^') , + oH, 
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and this gives the unique solution 

(4.19) Y = nB - 1 [(ln|l-/ 2 |/o / )'] + o(n) as n -> 0. 

One can see that 

(ln|l-/^0 / GL;(R + ) > 

so (I4.19p makes sense. Thus, (I4.12p and (14. 19|) actually determine the behaviour of the 
n-branch of similarity profiles of Riemann's problem for sufficiently small n > 0. 



4.4. On global continuation of n-branches: open problem. Global continuation 
of branches of nonlinear eigenfunctions of (14.11) from the branching point is often an 
intriguing open problem. Global bifurcation results concerning continuous branches of 
solutions originated at n = are already given in Krasnosel'skii (the first Russian edition 
was published in 1956), [Ml p. 196]. Concerning further results and extensions, see 
references in [191 Ch. 10] (especially, see [19, p. 401] for typical global continuation of 
bifurcation branches), and also [35l § 56.4]. 

In general, for the integral equation (I4.13P with compact Hammerstein operators in 
weighted L 2 -spaces, it is known since Rabinowitz's study (1971) that branches are infin- 
itely extensible and can end up at further bifurcation points; see [TjJl § 29] for further 
information stated in the framework of bifurcation analysis and [35, § 56.4]. The above 
n-extension of the branch is possible at any n = no > provided that the linearized 
operator has proper spectral properties. For oscillatory profiles f(y) about / = 1, with a 
complicated infinite set of intersection points, this is difficult to check in general but, in 
principle, can be done in specific weighted spaces, provided assuming that the oscillatory 
structure near interfaces is known in detail. Nevertheless, for the non-variational eigen- 
value problem (14.1 j) with non-divergent and non-monotone operators, a rigorous treatment 
of global behaviour of branches is very difficult and remains open. 

Since we have obtained the unique continuous branch originated at n = 0, it will never 
come back to the origin, so this suggests that it can be extended up to n = 1 and even 
further until homoclinic and other bifurcation points will destroy its necessary quality. 
Nevertheless, we do not know principles of such a global continuation, so we end up our 
branching analysis as follows: 

Conjecture 0Jl. The n-branch of nonlinear eigenfunctions of fl^.ip that is originated at 
n = from /o exists for all n G [0, 1] and does not have turning {saddle-node) points. 

In Figure [7J we show the deformation with n, with the step An = 0.1, of the CP 
similarity profile fo(y) satisfying (14.11) . It is seen that all the profiles have quite similar 
shapes for n = (the dotted line) and n — 1 including also the case of the negative 
exponents 

n = -0.1, -0.2, -0.3, -0.4, -0.5, 

for which fo(y) does not have finite interface but remains equally oscillatory as y — > +oo. 
Figure [8] shows the enlarged oscillations about the constant equilibrium 1 of the profiles 
from Figure [7J 
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Similarity profiles of the CP for n e [-0.5,1] 
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FIGURE 7. Continuous deformation in n of similarity profiles satisfying (|4.ip ; 



Similarity profiles of the CP for n e [-0.5,1 ], enlarged oscillations 
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FIGURE 8. Enlarged oscillations of similarity profiles from Figure n £ [— ^, 1]. 



5. RlEMANN PROBLEM |: SIMILARITY SOLUTIONS 

We now consider the RP-| fTTTTTI) . ffl~19D . (QUI) or (OB . Key aspects of the similarity 
analysis (both theoretical and numerical) are very similar. We stress the attention to some 
distinctive features. 

In Figure [9j we show the unique similarity profile fo(y) of the RP-| in the Cauchy 
setting. Notice that at the equilibrium level / = 1 the similarity profile is clearly less 
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CP: m=2, n=1 , the unique profile f (y) 




FIGURE 9. The CP solution of the RP-± (fTTTl) . (091) . (fL2TD for n = 1. 
Shooting FBP-1/2: m=2, n=1, first profiles f (y) and f (y) 
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FIGURE 10. Two first FBP solutions of the RP-± (fTTTD . (TTTT9T) . (TOUD for n = 1. 

oscillatory than at {/ = 0}, where propagation is approximately governed by the linear 
bi-harmonic PDE (fT~24f) . 

Figure [10] shows in the enlarged scale first two FBP profiles for this RP— ~ that are 
rather close to the CP one. The interfaces are 

yoi ~ 2.07 and yo2 ~ 3.36. 

The key aspects of homotopic connections and branching at n = for the RP-| remain 
the same as for the RP-1. The only change is that, convo luting with Heaviside data (11.151) 
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yields for n = the following similarity profile: 

+00 

(5.1) f (y)= j F(z)dz (n = 0). 

-y 

6. Interface equations for the TFE-4 with unstable linear diffusion 
We consider the TFE ( ll.33p . and, studying finite propagation close to u — 1, we set 

1 — u = v 

to obtain, up to small perturbations, the following unstable TFE: 

(6.1) v t = -(\v\v xxx ) x - v xx . 

As usual, we take the simplest travelling wave (TW) solutions 

(6.2) u(x,t) = f(y), y = x-\t =► -A/' = -{\f\f")' - /"• 

Hence, on integration once assuming the zero flux condition at the interface at yo = 0, 

(6.3) Xf=\f\f" + f for y>0, /(0) = 0. 

6.1. The FBP coincides with the CP. In the first approximation, we neglect the 
non-stationary term on the left-hand side of (I6.3P and consider the equation 

(6.4) |/|/"' + /' = for y>0, /(0) = /'(0) = 0. 

If /'(0) = 0, then the flux vanishes due to this ODE. Taking the expansion f(y) = y 2 p(y) 
and substituting into (16 .4p we get 2p' = — 1 + ... . This gives the following first term 
expansion: 

(6.5) f(y) = y 2 \\ny\ + .... 

Consider next the full equation (16.31) . Dividing by |/| and integrating yields 

(6.6) X y = f" + \n\f\+C (C EM). 

Proposition 6.1. (i) Equation (Iff. 61) does not admit oscillatory solutions near the inter- 
face at y = 0; and 

(ii) (1 6'. 6]) does not have a decaying solution with /(+00) = 0. 

In fact, (ii) shows that (16.61) does not admit infinite propagation (for TWs of course). 

Proof, (i) Indeed, assuming that f(y) — > as y — * 0, we have ln|/(y)| —1, so that 
f"{y) > and f(y) is strictly convex in a neighbourhood, (ii) is straightforward. □ 

This means that the asymptotics (I6.5P is already the actual maximal regularity be- 
haviour corresponding to both the FBP and CP, which thus coincide for such a PDE. 

The dynamic equation for the interface at x — s(t) follows from (16. 6p on differentiation, 
so we obtain the following equation with a third-order interface operator S: 

(6.7) s '(t) = \ = (f" + ln\f\)' = S[u]=v xxx + f at x = s(t). 
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A rigorous justification of such equations assumes using the von Mises variable X(v,t) 
locally near sufficiently smooth and monotone interface, 

X(v(x, t),t) = x, 

though it leads to a number of technical difficulties; see 



7. Interface equations for the unstable nonlinear diffusion 

We now treat the interface propagation for the equation (ll.35p . which, for v = 1—u ~ 0, 
reduces to the standard TFE-4 with the unstable diffusion 

(7.1) v t = -(\v \v xxx ) x - {\v\ m v x ) x . 
The TWs as in (16. 2p . on integration, yield the ODE 

(7.2) A/ = i/ir+i/r/ ; - 

7.1. The FBP. Case I. The local analysis close to the interface at y — yo is pretty 
standard. Namely, we have got the pure TFE expansion like (12.21) [TOl 



(7.3) f{y) = Cy 2 + C 1 y 3 + ..., where C\ = |, C > 0, 
provided that 

m > |. 

This defines the standard interface equation 

(7.4) s'(t) = A = 6Ci = v xxx at x = s(t). 

Case II. For the critical exponent m = |, all three terms in {\7.2h are involved in the 
expansion yielding in (17.31) 

Ul _ 6 _ 3 ° ' 

Therefore, the interface equation for TWs contains two operators 

(7.5) s'(t) = A = S[u] = v xxx + 2(| v xx ) m at x = s(t). 

Case HI. Finally, for m e (0, \), the FBP expansion is not C 3 -smooth. For instance, in 
the first such range m G (7,3)? the expansion is 

(7.6) f(y) = Cy* + C iy ^ + ly3 + ..., where C x = - 2m(1+ 4 (1+2m) C™, 

and, in the second term, the exponent satisfies 2 < 2 + 2m < 3. Therefore, / ^ C 3 ([0, 1]), 
so the equality A = /"'(0) makes no sense and the interface equation becomes rather tricky. 
The simplest way to resolve (17.61) with respect to the speed A = s'(t) is to decompose the 
TW solutions as follows: 

v(x, t) = v(x, t) + v(x, t), 

where v(x, t) is C 3 -smooth at the interface and v(x, t) is not. Then the interface equation 
in this non-regular case takes the standard (similar to (17.41) ) form 

s'(t) = {v) xxx at x = s{t). 
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An alternative, higher-order version can be obtained by observing from (17.21) that, for 
solutions (17. 6p . 

\Cy 2 + ... = ff'" + f m f (C=i/"(0)). 
Differentiating two times yields the following interface equation: 

s'(t) = S[u} = ^ {vv xxx + v m v x ) xx at x = s(t). 

The case m = \ is critical, where the second term on the right-hand side in (I7.2p adds an 
extra operator into the interface equation that now reads (the above differential form is 
also available) 

s'(t) = {v) xxx - l{\v xx ) 2m 1 at x = s(t). 
A rigorous proof of such types of interface equations is a difficult open problem. 

7.2. The CP. Case I. For m > ±, the propagating (A < 0) TW solutions are oscillatory 
near interfaces and have the behaviour as in ( 13. ip . i.e., 

(7.7) f(y) = yMs), s = lny, 

where (p(s) is a periodic solution of the ODE ( 13.21) with Ao i— > —A. 

Case II. In the critical case m — |, we still have the expansion (17.71) . where the oscillatory 
component solves a different ODE, 

(7.8) P 3 ((p) = A sign y?- (<pf + 3(p)\(p\~* . 

Existence-uniqueness of a stable periodic orbit becomes harder and remains an open 
problem. Figure [TTTa) shows the unique stable periodic solution of (17. 8ft . In Figure (b), 
for comparison, we show the periodic behaviour in the ODE for the stable counterpart of 
the TFE with m = |, 

1 2 

(7.9) v t = -(\v\v xxx ) x + (\v\3v x ) x =>> P 3 ((p) = Asigny? + (cp' + 3cp)\(p\~3. 
Case HI. Finally, in the range 

< m < min { ^ , n| 

two terms on the right-hand side of (I7.2p are leading that yields the same quadratic 
expansion (17. 6p as in the FBP. It seems that, in this parameter range (and also for m = 
in Section [6]), the CP and FBP settings coincide. 

In the CP, in view of the oscillatory nature of solutions in most cases, even formal 
derivation of interface equations becomes difficult, to say nothing about their rigorous 
justification. 

8. On higher-order TFEs: the RP-1 for m = 3 

Similarity solutions persist for higher-order TFEs (ll.32p . The matching becomes more- 
parametric and the topology of local bundles gets rather delicate. Numerics also get 
essentially more involved. 
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(a) unstable diffusion, m 



(b) stable diffusion, m 



FIGURE 1 1 . Stable periodic changing sign behaviour of the oscillatory compo- 
nent for the ODE (jEHJ (a), and (JEIl (b); A = — 1. 



8.1. On the CP profile and oscillations. As an example, in Figure [T21 we present the 
CP profile satisfying (11.321) for the case m = 3 and n = 1, i.e., for the sixth-order TFE 
(TFE-6). Instead of (13. ip . the oscillatory behaviour about f(y) = 1 is smoother, 

(8.1) f(y) = 1 - (y - y)V(s) + - , s = \n(y - y), 

where ip(s) solves a fifth-order ODE 



(8.2) P 5 (<p) = V (5 > + 15</ 4) + 85</' + 255</ + 27V + 12(fy> = A sign <p, 

where Ao = \ yo- See details on derivation in [25, § 12.2]. Figure [T51 shows the unstable 
periodic solution of (18.21) for Ao = 1, which by (18. ip describes the oscillatory behaviour 
near the interface at 



13) 



2/o,cp ~ 13. 



8.2. Shooting FBP profiles. The results of a standard shooting of the first FBP profiles 
are presented in Figure dH Similar to shooting in Figure [U here we have chosen the same 
strategy for higher-order equations: for m = 3, we fix three free-boundary conditions 
from (I1.33P excluding the second one, so changing the interface location yo as the free 
parameter, we want to satisfy the zero contact angle condition (I2.54p . As a consequence, 
we see in Figure [T3] a consistent illustration of the limit FBP-CP process indicated earlier 
in the schematic Figure [6j Namely, we obtain four first FBP profiles, fi(y), f2(y), f3(y), 
and fi(y), with interfaces positions at 

(8.4) ?/oi ~ 3.05, y 02 ~ 4.85, y 3 ~ 6.55, and y 04 ~ 8.1. 

These are smaller than the Cauchy one (18. 3p . which is expected to be the limit of all other 
FBP interfaces {yok} as k — > oo. 

As usual, only the first profile fi(y) satisfies the physical range condition (11.231) . Other 
profiles are oscillatory about the constant equilibrium / = 1, and as it is seen in Figure 
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Similarity profile of the CP for n=1 , m=3 



Similarity profile of the CP for n=1 , m=3; enlarged oscillations 




10 12 14 



(a) CP profile 



(b) oscillations, enlarged 



FIGURE 12. The unique CP profile f(y) solving (fQ2l) . (fOTD for m = 3, n = 1; 
the profile (a), and enlarged oscillations about 1 (b). 



periodic orbit: m=3 



-e- 




1.5 2 2.5 3 3.5 4 



FIGURE 13. The trace of the unstable periodic orbit of ()8.2j) . with Ao = 1 (yo = 6). 

[T4l rather thoroughly mimic the oscillations of the linear solution (the boldfaced dotted 
line) defined by the same formula (ll.28p . This figure confirms that Conjecture [31 2 applies 
also to the TFE-6 (and seems to any of 2mth-order ones (11.301) ). 

For comparison, in Figure [151 we present first three FBP profiles for the linear case 
n = 0, i.e., / solves the FBP 

/ (6) + |yf = for y>0, 
/(0) = /"(0) = /( 4 )(0) = 0, f(y ) = 1, f(y ) = f"(y ) = f^(y ) = 0. 

The interface positions are larger than that in (18. 4| for n = 1, 



yoi ~ 3.29, ?/02 ~ 6.57, and y 03 ~ 9.72. 
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Shooting FBP profiles for n=1 , m=3 
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FIGURE 14. Shooting first four FBP profiles satisfying (fL32l) . (Q3l) for m = 3, 
n = 1. 

Observe that the FBP profiles are close to corresponding humps of the CP profile denoted 
by the boldface dotted line. The general geometry of the mutual location of FBP and the 
CP profiles for n = is quite similar to the "nonlinear" one in Figure dH for n = 1 (the 
figures are topologically equivalent). 

8.3. Branching at n = 0. For n = 0, (jl.30p for m = 3 becomes the tri-harmonic equation 

(8.6) u t = u xxxxxx in Kx R + . 
Its fundamental solution is given by 

(8.7) b(x,t)=t-$F(y), y = ^k, 
where F is the unique symmetric solution of the problem 

(8.8) BF = F^ + lF'y + lF = in R, J F = 1. 

Extension of the similarity profiles of the ODE (jl.32p from the branching point n = is 
the same as that for the TFE-4 in Section 14.11 Spectral properties of the corresponding 
linear operators B in (18.81) and the adjoint one 

B* = Dl-\ yD y 

in Lp(R) and LL(R) respectively can be found in [20J. 

In Figure [T6| we show branching of the similarity profiles at n = (the dotted line), 
where, with the step Ap = 0.2, we cover the range p G [—0.8,1]. All the profiles are 
oscillatory, but for n < the interface is situated at yo = +oo. 
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Shooting FBP profiles for n=0, m=3 

1.15 n 1 ; 1 1 1 1 




y 

FIGURE 15. Shooting first three FBP profiles satisfying the linear FBP problem (18. 5p . 



Similarity profiles of the CP for n e [-0.8,1]; m=3 




5 10 15 20 



y 

FIGURE 16. Similarity profiles of the CP satisfying (jl.32p for m = 3 appearing 
at the branching point n = (the dotted line). 
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